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We show from covariant transport theory that, for a massless ideal gas equation of state, even a 
small shear viscosity to entropy density ratio r\ m s/(4ir) generates significant 15 — 30% dissipative 
corrections to elliptic flow for conditions expected in mid-peripheral (6 = 8 fm) Au + Au collisions 
at ^/~s17n ~ 200 GeV at RHIC. 
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PACS numbers: 12.38.Mh, 25.75.-q, 25.75.Ld 

Introduction. The goal of the heavy-ion physics is to 
study the properties of nuclear matter at extreme en- 
ergy densities and temperatures. Spectacular features 
of the data from the Relativistic Heavy Ion Collider 
(RHIC) for gold-gold reactions at center-of-mass energies 
y/sNN ~ 100 — 200 GeV lead to the suggestion that the 
hot and dense matter created in the collision exhibits per- 
fect fluid behavior. (A perfect fluid has vanishing shear 
and bulk viscosities and heat conductivity.) The cor- 
nerstone of this conclusion was the success of ideal (Eu- 
lcr) fluid dynamics in explaining the large second Fourier 
moment of the azimuthal momentum distribution for 
fixed vt, the so-called "elliptic flow" V2(pt) = (cos2^) PT 
[E S S SI [H ; observed in non-central collisions. 

General considerations based on quantum mechanics, 
on the other hand, indicate a nonzero lower bound for 
shear viscosity. For example Q, combining the time- 
energy uncertainty principle with kinetic theory gives 
rj/s > 1/15, where s is the entropy density. This result 
is also supported by calculations for Af = 4 supersym- 
metric Yang-Mills theory via the gauge-theory - grav- 
ity AdS-CFT duality conjecture. Pioneering results pre- 
dicted rj/s > 1/(4tt) 0, while recent calculations for an 
extended class of theories find a somewhat lower bound 
rj/s > 16/25 x l/(47r) [8(. It is an open question whether 
either of these limits applies to QCD. But the possibility 
is intriguing because even a small rj/s ~ 0(1) /(4tt) has 
significant dynamical effects in heavy-ion collisions. 

Ideal hydrodynamics assumes local equilibrium 
throughout the evolution. For nonzero transport 
coefficients, on the other hand, the system departs 
from local equilibrium, leading to dissipative correc- 
tions. If the system stays sufficiently close to local 
equilibrium, dissipation can be investigated via causal 
dissipative hydrodynamics, for example, Israel-Stewart 
theory [rJ, EH- Solution techniques for viscous 
hydrodynamic equations, in the minimally required 
2+1 dimensions necessary for elliptic flow studies, 
have been recently developed and a pplieddl, m, m. 
However, Israel- Stewart theory comes from a truncation 
procedure 0, [n| with uncontrolled errors (it lacks a 
small expansion parameter), and therefore its region of 
validity is not known. Moreover, the solutions are causal 



only in a region of hydrodynamic parameters, and their 
stability is not guaranteed[16|. 

Here we utilize instead the full y c ausal and stable co- 
variant parton transport theory [171, lis llE 120, [21| for 
which covariant algorithms have been available in full 
3+ ID for quite some time[22l. liij . I24I ]. A forerunner of 
this analysis 2l| considered large, constant elastic 2^2 
gluon-gluon cross sections a ~ 45 mb, and found a sig- 
nificant 20 — 30% reduction of elliptic flow due to dis- 
sipation. Such dynamics gives an increasing rj/s with 
time (25)]. In contrast, here we study dissipation for a 
"minimal" -q/s « l/(47r) that is constant in time. Pre- 
liminary results have been reported already at [26j, and 
are confirmed here to be accurate. 

Covariant transport theory near the hydrodynamic 
limit. We consider here, as in Refs. [H, El, EI 0,113, 
the simplest but nonlinear form of Lorentz-covariant 
Boltzmann transport theory in which the on-shell phase 
space density f(x,p), evolves with an elastic 2^2 rate 
as 

Pid^h = S(x,fi) + \JJJ (hh - hh) 

234 

X S 4 (pi+P 2 -P3-Pi) ■ (1) 

where the integrals are shorthands for J. = J d 3 pi/(2Ei). 
For dilute systems, / would be the phasespace distribu- 
tion of quasi-particles, while the transition probability 
W = s(s — Am 2 )da / dt would be given by the scattering 
matrix element. Our interest here, on the other hand, is 
to study the theory near its hydrodynamic limit. 

It is well known (Boltzmann's H -theorem) that ([T]) 
drives the system towards a fixed point, global equi- 
librium. In the hydrodynamic limit (W — > 00), the 
transport solutions approach local equilibrium f(x,p) — 
gexp[(fi(x) — pvu" (x))/T(x)]/ (2ir) 3 . A systematic ex- 
pansion in small gradients around equilibrium via the 
Chapman-Enskog procedure [l5| gives the viscous hydro- 
dynamic equations by Navier and Stokes. However, this 
approximate theory is severel y a causal. A causal for- 
mulation proposed by Mueller 27| and later generalized 
by Israel and Stuart (IS) [§] retains certain second-order 
derivative terms, resulting in relaxation equations. IS 
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theory can also be recovered from transport via the 14- 
moment expansion of Grad [tj [TBj]. The transport co- 
efficients, and the microscopic relaxation times for the 
dissipative fluxes in IS, are all given by the differential 
cross section da/dt. 

The key observation here is that one can use trans- 
port theory to solve causal viscous hydrodynamics pro- 
vided one dials in the equation of state (EOS) and trans- 
port coefficients of interest. In this case, the "particles" 
and the specific "interaction" in the transport have no 
physical significance - they are only mathematical tools 
to reproduce the desired dynamical equations. Here we 
consider an ultrarelativistic gluon gas with e = 3p, appli- 
cable to the high-temperature plasma in the early stages 
at RHIC. In this cased Q> V ~ 4T/(5cr tr ) and the 
shear stress relaxation time is r ff = 6At r /5, where o~t r and 
At,- = l/(no~t r ) are the transport cross section and trans- 
port mean free path, respectively. Note, for an isotropic 
cross section cr tr = 2cr to t/3. 

Elliptic flow and rj/s. For the above conditions, 



JL 

An 



TXt 



T 



5na± 



(2) 



Assuming the system stays close to local equilibrium, 
during the initial longitudinal (Bjorken) expansion stage 
of the heavy-ion collision the density and temperature 
evolve as n ~ 1/r, T(t) ~ r -1 / 3 where r = \Jt 2 — z 2 is 
the longitudinal (Bjorken) proper time. For a constant 
cross section, rj/s oc t 2 / 3 then increases with time. 

One might therefore think that Ref. [2l| with constant 
a ~ 45 mb overestimated dissipative effects in Au + Au 
at RHIC energies. However, even though rj/s grew in 
that calculation, its initial value was really small (cf. 
Fig. [TJ. For the longitudinally boost invariant scenario 
assumed there, Xt r — t /{o~tr dN /dndx 2 -, ) , and with the 
parameters of that calculation [31| Af r (ro = 0.1/to) ps 
7.1 x 1CT 3 fm(!), i.e., rj/s ~ l/(607r) at the very center 
of the collision. The ratio is way below the conjectured 
bounds even for the average density that is 2 — 3 times 
smaller than the maximum. 

We cross-check this important finding with the trans- 
port opacity [2oT] 



X = {ncoii} = 

0~tot 



dz 



A tr (x + zn, t = r + z) 



(3) 



which is the number of collisions per particle weighted 
by the transport cross section and averaged over initial 
coordinates and directions, x is dominated by the early 
and densest longitudinal expansion stage, during which 
Xtr fx r, and thus 



X 



dz 



(Xtr (to)) Jo z + t 



to L 
(Xtr (to)) Tq 



(4) 



With x ~ 21 from the calculation and an estimate L ~ 
3 — 4 fm for the size of the reaction zone, on average 



(A tr (ro = 0.1/to)) ps 1.5 x 1(T 2 fm. This is about 2.5 
times larger than the value estimated for the collision 
center, as expected. 
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FIG. 1: Schematic evolution of n/s for the calculation by 
Ref. [2l( | that considered a constant cross section. 

From ^ we then find that for the situation in Ref. [2l[ , 
rj/s evolves with time schematically as shown in Fig. [TJ 
During the first few fermis of the evolution relevant for 
the buildup of elliptic flow [2l[, the system stays 
closer to equilibrium than would be allowed by a "min- 
imal" viscosity because the transport mean free path 
(or equivalently, the scattering rate) is not limited by 
any quantum bound. The situation only changes after 
t ~ 1 — 3 fm, when the transport mean free path be- 
comes large enough to reach rj/s ps l/(47r). 

To study elliptic flow for constant rj/s, we therefore 
start with an initially modest transport cross section 
and increase it with time as atrir) — <7o[t/(0.1/to)] 2 ' 3 . 
Such growth is already encoded in perturbative QCD, 
provided we ignore the running of the coupling: a tr ~ 
(187ra 2 /s)ln(s/M!>) - (™ 2 JT 2 ) m(18/ 5 2 ), where ji D - 
gT is the Debye mass. The initial conditions for Au + Au 
at ^snn = 200 GeV at impact parameter b = 8 fm 
(see footnote [31]) and the numerical solution technique 



MPC[23j are the same as in Ref. [21(. For numerical 
convenience we uniformly set o~o ~ 2.7 mb, which en- 
sures that on average rj/s ps 1/(4tt) in the system. In 
the center of the collision zone, ([2]) gives a lower value 
rj/s « 0.4/ (47r), but that is compensated by the increase 
of rj/s with decreasing density as we go outward. With 
the average density (n) ~ n max /2.5 estimated earlier 
from f4|, on average Airrj/s pa 1. A cross-check with 
the transport opacity x « 16 obtained from the growing- 
cross-section calculation er tr (r) cx r 2 / 3 gives 
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i.e., X tr (T =0Afm) - 0.09-0.11 fm and inrj/s - 0.8-1. 

The above choice of cro implies ~ 5 — 10 times higher 
two-body rates than perturbative QCD estimates. For 
our purposes this is not a problem, the rates are simply 
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adjusted to reproduce the desired rj/s. With the inclusion 
of radiative 3 <-*■ 2 processes [29j], even perturbative rates 
could generate a small rj/s ~ 0.1. 

With 2 — > 2 scattering, particle number is conserved 
and thus s = n(4 — /x/T) where jx is the chemical po- 
tential. Therefore, @ acquires a small relative correc- 
tion As/s(M=0) = (l/4)ln(n eg /n) = (1/4) h%T 3 /(7r 2 n)] 
logarithmic in density. For our gluon gas initial condi- 
tions (g = 16), it is about ~ (—6)% at the center of 
the collision, while ~ 20% for a low n — n max /3. Dur- 
ing the longitudinal expansion stage, the correction stays 
roughly constant because the dilution is largely compen- 
sated by cooling, nT 3 w const. Dissipation of course 
still generates entropy because the temperature drops 
slightly slower^ than T oc r -1 / 3 , but that effect cannot 
be very large for the system to stay near equilibrium. 
With Airn/s ~ 0.8 — 1 initially, we have a cushion for 
entropy production during later evolution. Therefore, 
we conclude that the averge rj/s is set to the desired 
rj/s = 1/(4-71") within about 20% in the calculation. 

Figure [2] shows differential elliptic flow V2(pt) results 
for To = 0.1 fm. Even for a "minimal" rj/s = l/(4w) 
(filled squares), dissipation reduces elliptic flow at mod- 
erate pr ~ 2 — 3 GeV by about 25% relative to the ideal 
hydrodynamic limit (solid line). The relative change in- 
creases with decreasing pt, and therefore dissipation also 
flattens the slope at low vt [3^j . For comparison, we also 
show the result (open squares) for a constant cross sec- 
tion a tr = Co, he., growing rj/s w (t/to) 2 / 3 /(47t), which 
of course generates much smaller elliptic flow. 



ideal hydro 
i - M 




■q/s « 1/(4tt) 



□ □ □ □ 



1 2 
Pt [GeV] 



FIG. 2: Differential elliptic flow V2{pt) as a function of pr in 
Au + Au at ^/s^F = 200 GeV and b = 8 f m at RHIC, from 
ideal hydrodynamics (solid curve) using the codes in @, 0] 
and covariant transport (squares) using the MPC algorithm 
[23l ]. An initial (thermalization) time of To — 0.1 fm/c was 
assumed. Transport results for a constant cross section (open 
squares) and for n/s « l/(47r) (filled squares) are shown, 
while the hydrodynamic curve is from [2lT |. 

Figure [3] compares differential elliptic flow 1^2 (pt) from 
the transport (squares) to the ideal hydrodynamic limit 
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FIG. 3: Same as Fig.[2]but for an initial (thermalization) time 
ro = 0.6 fm/c. 



(solid line) for a later thermalization time ro = 0.6 fm. 
In this case the initial temperature was adjusted to 
T w 0.385 GeV to account for cooling T ~ t^ 1 / 3 be- 
cause of pdV work during longitudinal Bjorken expan- 
sion. With the rescaled temperature, results for both 
ideal hydrodynamics and transport with constant cross 
section (open squares) are essentially independent of ro, 
as long as tq /R <C 1, as was also found in Ref. [2l| (cf. 
Fig. [2]). However, with rj/s « 1/(4tt), the dissipative re- 
duction of elliptic flow (solid squares) relative to ideal 
hydrodynamics is more modest, ~ 15%, for the larger 
to = 0.6 fm. This is similar in magnitude to recent re- 
sults from 2+ ID dissipative hydrodynamics 12j|l4|. 

At fixed 77/s, dissipative effects are weaker for larger 
To because the initial value of the transport cross section 
is larger in that case. An equivalent explanation is that 
though rj/s is the same, initial velocity gradients in the 



system d t 



1/r are smaller for larger tq. Indeed, the 



Navier-Stokes correction to the stress tensor [rl flFijj 



\A^d a u a ) (6) 



(A^ = g» v - u^u v , W = A^ v d v , is the flow veloc- 
ity), in case of a longitudinal boost-invariant expansion, 
implies viscous corrections to the transverse and longitu- 
dinal pressure Apx = 2?y/(3r), Api = — 4?//(3t). There- 
fore, relative pressure corrections 



Ap r/ Ts Amj 1 
p s pTr s ttTt 



(7) 



if rj/s — const; 



decrease with time Ap/p cx r~ 2//3 , 
whereas Ap/p oc r , if o tr — const. 

Based on the Navier-Stokes estimate <[TJ) , it is not a sur- 
prise that dissipative corrections are important for condi- 
tions expected at RHIC. Ap/p ~ 20% for r T = 0.6 fm 
x 0.385 GeV and is almost 100% for t T = 0.1 fm 
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x 0.7 GeV. In the former case, the correction is mod- 
est, and viscous hydrodynamics is likely applicable. In 
the latter case, however, a hydrodynamic approach seems 
questionable. It would be interesting to test this antici- 
pated break-down of hydrodynamics against solutions of 
Israel- Stewart theory. We note that dissipation is ex- 
pected to be relevant not only at RHIC but at the LHC 
as welllim. 

Conclusions. We utilized covariant transport theory to 
study the effect of a "minimal" shear viscosity n = s/(4V) 
on differential elliptic flow V2(pt) in Au + Au collisions 
at ^/sjvat ~ 200 GeV at RHIC. The key ingredient is 
a transport cross section cr tr oc t 2 / 3 that grows with 
time. We find significant reduction of elliptic flow rel- 
ative to the ideal hydrodynamic limit, ~ 25% reduction 
for a thermalization time To =0.1 fm, while ~ 15% for 
t = 0.6 fm. This indicates that even such a small shear 
viscosity cannot be ignored at RHIC, and thus the LHC 
as well, because gradients are large. 

We note that this study set n/s w l/(47r) only within 
20% and in an average sense, for numerical convenience. 
The evolution of the local density and temperature 
was approximated with analytic results for longitudinal 
Bjorken expansion. More accurate results could be ob- 
tained, in principle, via a transport cross section <7t r (n, r) 
that depends explicitly on the local density. However, 
that is much more expensive numerically with the solu- 
tion technique (cascade algorithm) utilized here. 

In principle, radiative processes such as gg <-> ggg can 
also be included [29|. However, near the hydrodynamic 
limit we do not expect large corrections to our results 
because the dynamics is determined solely by the equa- 
tion of state and the viscosity r)/s = 1/(4tt). The main 
difference is that radiative processes allow for change of 
particle number, which should be only a modest refine- 
ment for the initial conditions considered in this study. 

Finally, we emphasize that a simple ideal gas equation 
of state e = 3p has been considered, which also implies 
vanishing bulk viscosity. It would be important to repeat 
this study with an equation of state that is more realistic 
for quark gluon matter at moderate T < 300 MeV, and 
to investigate effects of bulk viscosity which is expected 
to rise sharply in the vicinity of T ~ 200 MeV[3C§. 
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